Transverse fluctuations of grafted polymers 
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Abstract 

We study the statistical mechanics of grafted polymers of arbitrary stiffness in a two-dimensional 
embedding space with Monte Carlo simulations. The probability distribution function of the free 
end is found to be highly anisotropic and non-Gaussian for typical semiflexible polymers. The 
reduced distribution in the transverse direction, a Gaussian in the stiff and flexible limits, shows 
a double peak structure at intermediate stiffnesses. We also explore the response to a transverse 
force applied at the polymer free end. We identify F-Actin as an ideal benchmark for the effects 
discussed. 

PACS numbers: 87.15.Ya, 87.15.La, 87.16.Ka, 36.20.Ey 
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Healthy cells require an efficient and complex transport network to carry out the over- 
whelming number of tasks that are needed to accomplish their function. This network, also 
known as the cyto skeleton, is formed primarily by filaments (actin ffiaments, microtubules 
and intermediate filaments), linked together by a large collection of accessory proteins Q|. 
A complete description of the structural and mechanical properties of these filaments is 
therefore essential in order to unveil the mechanical properties of the entire cell. Advances 
in the field have been significantly promoted by a unique set of optical and mechanical tech- 
niques which allow to visualize and manipulate single cytoskeletal filaments 3, 0, Q| and 
DNA P|. Fluorescence videomicroscopy [^| and nano- manipulation [3] can be conveniently 
used to obtain quantities as the distribution function of the end-to-end distance or the 
mechanical response to an external force in great detail and at the single molecule level. 
These quantities are amenable to a direct comparison with theoretical models. 

The main material parameter in the description of a polymer filament is its persistence 
length, £p. It is defined as the typical length over which correlations of the tangent vectors 
of the filament contour decay. Polymers are considered to be flexible when their persistence 
length is small compared to their total length L, or t := L/£p > 10. In this limit, they can 
be well described by the minimal model of the Gaussian Chain Q]. Polymers of biological 
importance, e.g. F-actin, are often semifiexible, meaning that their persistence length is 
comparable to their total length. While flexible polymers are dominated by entropic effects, 
the statistical mechanics of semifiexible polymers is strongly affected by their bending energy 
and the close vicinity of the classical Euler instability for buckling a rigid beam j^. 

The distribution function P(R) of the end-to-end vector R, a simple Gaussian for a 

n 

flexible polymer, is peaked towards full stretching and is completely non-Gaussian [10||. 
The mechanical response of a semifiexible polymer is highly anisotropic, depending on the 
direction in which the force is applied ll|. These findings result in bulk properties of 
solutions and networks that are completely different from the isotropic elasticity of flexible 
polymer solutions 0, Q|. In addition, the inextensibility constraint becomes crucial in 
determining the approach to full stretching upon the application of a force /, as reported 
by Marko and Siggia for double-stranded DNA. 

Here we investigate the mechanical and statistical properties of a single chain grafted 
at one end, a problem of direct relevance for force generation in cellular systems. The 
other end is either free, or subject to a constant transverse force, whose magnitude extends 
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into the non-linear regime. We restrict ourselves to a two-dimensional embedding space, 
since in most experiments, fluctuations in one direction are severely restricted, or cannot be 
observed. The generalization to a three-dimensional space is straightforward and will be 
reported elsewhere lljl. 

n 

We refer to the Wormlike Chain Model (WLC) introduced by Kratky and Porod |15| . 
In this framework, a polymer conformation is represented by a succession of segments ti, 
whose direction is tangent to the polymer contour at the zth segment. Since the polymer 
is assumed to be inextensible, all segments have a prescribed length a = L/N. The 
Hamiltonian is given by: 

N-l N 
i=l i=l 

where e is the energy associated to each bond and / is a force eventually applied to the 
second end. It is also possible to deflne a continuum limit for a ^ 0, ^ oo, with Na = L 
and e = ea^jN held flxed. The Hamiltonian in Eq. ^ is then equivalent to the following 



functional 
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where k = eL and is the tangent vector of the space curve r(s) parametrized in terms 
of the arc length s. The inextensibility of the fllament is imposed by the local constraint 
= 1. The continuous version of the wormlike chain has been successfully used to 
obtain various statistical quantities, as the tangent-tangent correlation function or moments 



16, Il8||. It has been recently used to obtain the 
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[l^ , and force-extension relations 0, [ll[ [3l 



of the end-to-end distance distribution 
radial distribution function 

We have developed a Monte Carlo simulation to investigate the behavior of a semiflexible 
polymer in the proximity of the limit t — > 1. The rationale behind this choice is the search 
for clear hallmarks of the onset of the "semiflexible" nature of a fllament. In this inter- 
mediate limit, analytical results are difficult to obtain: typical approximation schemes that 
build on either Gaussian chains or rigid rods are outside their validity range; hence, com- 
puter simulations become crucial. The first end of the filament is assumed to be clamped, 
i.e. the orientation of its tangent vector is held fixed along a direction, named the x-axis. 
The second end is left free to assume any possible orientation. The initial configuration has 
been randomly chosen in the proximity of the full stretching condition, thus ensuring a fast 
convergence to equilibrium. A new configuration is generated by changing the orientation of 



one segment and accepted according to the standard Metropolis algorithm and the discrete 
Hamiltonian, Eq. ^ Effects resulting from self-avoidance are not considered, but we notice 
that configurations where the chain folds back onto itself are strongly energetically sup- 
pressed for sufficiently stiff polymers. Results ceased to depend on the number of segments 
for = 50. On the order of 10^ Monte Carlo steps per segment were performed, and results 
were averaged over different runs, obtaining a perfect agreement between measured expec- 
tation values of the end-to-end distance {R^) and {R^) with known exact expressions. The 
radial distribution function was calculated and coincided with the analytic results in 3| 
within the accuracy thereby reported. 

Here we are interested in the probability distribution function P(x, y) of the free end 
in the plane determined by the direction of the clamped end (x-axis) and the transverse 
one (y-axis). This quantity is directly accessible to experiments allowing for a quantitative 
comparison with our predictions. We will also consider the reduced distribution functions 
P{x) and P{y), obtained by integrating P{x,y) over the variables y and x, respectively. 

It is important to notice that when both ends are free, the radial distribution function is 
rotationally invariant and is therefore only a function of the distance R between the ends. 
Clamping one end breaks rotational symmetry and leads to distinctly different longitudinal 
and transverse distribution functions, P{x) and P{y). Nonetheless, the broken rotational 
symmetry does not affect the total energy of the configuration. This implies, and is in fact 
confirmed by our simulations (data not shown), that the longitudinal distribution function 
P{x) coincides with the radial distribution function P{R) of the end-to-end distance, apart 
from a constant normalization factor. The characteristic feature of this function is a crossover 
from a universal Gaussian shape centered at the origin with a characteristic width determined 
by the radius of gyration, to yet another universal shape whose peak is shifted towards 
full stretching and whose width is determined by a new longitudinal length scale L\\ oc L'^/ip. 

This has to be contrasted with the transverse distribution function. Not surprisingly, 
given the intrinsic isotropy of flexible polymers, the distribution P{y) is a Gaussian and 
identical to P{x) for high values of t. In the stiff limit, P{y), at variance with P{x), 
is again a Gaussian centered at y = 0, whose width is now given by a new transverse 

length scale = ^2L^/3ip 01 

|2l|. Surprisingly, at intermediate values the probability 

distribution function is not a smooth interpolation between these two Gaussian limits but 
shows interesting and qualitatively new features. As t approaches the value 1 from above 
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FIG. 1: Distribution function for the projection of the free end along the transverse direction, P{y), 
obtained by Monte Carlo simulations. Lengths are measured in units of L. Errors are comparable 
to the point size in the insets, (a) Appearance of double peaks for t < 2.5. (b) Re-entrance from 
the double peaks to a flat distribution in the stiff limit t < 0.75. Insets show details of the crossover 
regions. 

(flexible side), the Gaussian peak is first smeared out into an intrinsically non-Gaussian 
flat distribution (see Fig. ^). At t = 2.8 (see inset), the distribution contains three local 
maxima, but as t is decreased, the central peak at y = loses weight to the two symmetric 
peaks off the x axis. The double-peak structure is most pronounced around t ~ 1.5, i.e. 
L ^ IMp. 

As the stiffness is increased, P{y) recovers its fiat structure, as shown in Fig. ^p. Notice 
also (inset of Fig. ^) that at t = 0.75 the two peaks start to compete with a growing peak 
centered at y = 0, such that one finds a triple maxima shape again. Although intrinsically 
non-Gaussian, this central peak will eventually tend to a Gaussian distribution in the stiff 
limit. The re-entrance from the double peak structure to a fiat distribution is a genuine 
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hallmark of semiflexibility. This effect cannot be explained by analytical calculations using 
a harmonic (or weakly bending rod) approximation, whose prediction for P{y) would be a 
Gaussian centered at j^]. Higher order cumulant expansions about a Gaussian distribution 
have also failed to provide a fast convergence to our P{y)- An entirely analytical solution 



can be provided by the eigenfunction approach described in j22| for persistent random walks, 
although the connection to our probability distributions would only be numerical. 

Finally, let us emphasize that the double-peak structure of P{y) does not indicate a 
bistability in the constant force ensemble. As shown below, linear response theory leads to 
positive force constants in this regime. What actually happens under the application of an 
external force is that the distribution function becomes asymmetric and weight is shifted 
from one peak to the other. In an experimental setting with a fixed transverse distance y 
and a correspondingly adjusting force, one would probe P{y) directly and be able to observe 
a kind of "bistability" . 

Further insight can be gained by the inspection of the joint distribution function P(x, y), 
represented with density plots in Fig. |21 In the stiff limit, y) should be confined to the 
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FIG. 2: Density plots obtained by Monte Carlo simulations: dense regions are colored in red, 
scarcely populated in blue on a color scale appropriately chosen to enhance the contrast, (a) t = 
2/3; (b) t=l]{c)t = 2; (d) t = 20. 



classical contour obtained by applying the elasticity equations to a rigid rod. This contour 
can be approximated by a parabola in the proximity of full stretching and is obtained through 
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elliptic functions for any deformation j9| . In Fig. |2K the classical contour coincides with the 
ridge of the probability distribution function. As we relax the stiffness, thermal fluctuations 
will make the tip of the filament explore the configuration space in the vicinity of the classical 
contour. Roughly speaking, transverse (bending) fluctuations enhance fluctuations along the 
classical contour and shift weight from the center to the upper and lower wings in Fig. l^tL-b. 
In contrast, longitudinal fluctuations widen the distribution function perpendicular to the 
classical contour. Since for a semiflexible polymer, the corresponding lengths L\\ and L±_ 
scale differently (transverse fluctuations are much "softer" than longitudinal ones), upon 
lowering the stiffness P(x, y) gains more weight in the wings rather than in the center. It is 
precisely this effect that gives rise to the double peak distribution, when P(x, y) is projected 
in the transverse direction (see Fig. |2b)- Eventually, in the flexible limit, where transverse 
and longitudinal fluctuations become comparable, P(x, y) is spread so as to cover almost all 
the available space (Fig.|2fc), before the isotropic Gaussian distribution is recovered (Fig.|21i). 

We have also explored the transverse response of semiflexible polymers by applying a 
constant force / in the transverse direction. The effect of a small applied force on the 
average end-to-end distance (or force extension relation) has been studied within linear- 

n 

response in In this work, we will consider the effect of an external transverse force of 
arbitrary magnitude on the average position {x) ^ and {y) ^ of the free end. 

In general, we expect {y) j to have the same parity of the applied force, and hence to be 
odd, while should not depend on the sign of the force and hence should be even. In 
the continuum limit, it is possible to write down the exact expressions for and {y) ^ and 
to show that the expected parities hold on very general grounds and that the response of 
the longitudinal extension to a transverse force is intrinsically nonlinear in the small force 
regime. Monte Carlo simulations confirm these predictions, as shown in Fig. El The response 
in the direction of the clamped end is even in / and it can be approximated by a parabola 
centered on the / = axis. The response in the transverse direction is odd in / and shows 
the same re-entrance phenomenon reported in |Q] for the linear response coefficient. 

Note that while in the case of a longitudinal force, the approach towards full stretching 
(or saturation) can be calculated within the weakly bending rod approximation, this is 
no longer true for transverse forces. The position of the free end can be calculated from 
classical elasticity theory P] and expressed by means of elliptic functions. Only in the high 
force regime or in the stiff limit, when fluctuations become unimportant, results from our 
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simulations coincide with classical elasticity theory. 
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FIG. 3: Response to a transverse force, obtained by Monte Carlo simulations. Forces are measured 
in units of ksT/L, lengths in units of L. Error bars are shown. (Above) Response in the clamping 
direction. (Below) Response in the transverse direction is odd with /: only part of the explored 
parameter region is shown for clarity. 

The effects hereby reported are amenable to a direct comparison with experiments regard- 
ing cytoskeletal filaments, or even DNA. For instance, optical systems might be used to get 
the X OT y projection of the radial distribution function for a particular class of semifiexible 
polymers. For F-Actin with ip ^ 16 fim Q], the double peak effect should be well visible for 
a range of lengths, 12 fim ^ L < A3 fim. In this parameter range the difference between the 
central relative minimum and the double peaks maxima results in 10% of the total length 
(see Fig. in the range 1 -h 4 fim that is well above the experimental precision of 0.05 fim 
reported by 0]. Hence F-Actin would provide an ideal benchmark for the effects we report. 
We emphasize that the double peak structure is a clear hallmark of semifiexibility and hence 
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it might be used to obtain a rough estimate of the persistence length of a particular polymer 
filament, as for instance the nanometer sized stalks of kinesins and myosins. 

In summary, we have presented evidence from extensive Monte Carlo simulations that the 
parameter region corresponding to semifiexible polymers is hallmarked by the appearance 
of a series of effects in the radial distribution function and in the response of the clamped 
polymer to an external transverse force. A semifiexible polymer shows a distinct anisotropy 
in the probability distribution function of the free end along the direction of the clamped 
end. At intermediate stiffness, L ^ ip, the distribution function shows a pronounced dou- 
ble peak structure in the transverse direction. Semifiexible polymers have been previously 
reported to be anisotropic objects, i.e. to respond in different ways to forces applied 
in the clamping or transverse direction. Here we have shown that even their response to 
a force along the transverse direction alone is intrinsically anisotropic, being linear in the 
transverse direction and non-linear along the direction of the clamped end in the small force 
regime. 
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T. Franosch and K. Kroy. This research has been supported by a Marie Curie Fellowship 
under contract no. HPMF-CT-2001-01432. 
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